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Abstract 

We investigate the deformation properties of atomic nuclei in a hadronic chiral SUfl avor (3) model 
approach. The parameters are fitted to hadron mass properties and adjustments for spherical 
finite nuclei have been performed. Using these parameters the deformation of a series of light and 
heavy nuclei are obtained in a two-dimensional self-consistent calculation. In addition a case of 
superdeformation in a heavy nucleus is studied. 
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I. INTRODUCTION 



Relativistic descriptions of strong interaction physics using hadronic degrees of freedom 
have again become more popular in recent years. One of the areas where hadronic models 
have been studied is the regime of highly excited hot and dense nuclear matter as explored 
experimentally using ultrarelativistic heavy-ion collisions. Towards the phase transition to a 
chirally restored phase a variety of models loosely connected to the chiral a model have been 
investigated. Extensions including vector meson fields (er — uj model) and the strangeness 
degree of freedom by incorporating hyperons and strange mesons have been discussed. In 
that context density and temperature dependent hadronic masses and the structure of the 
phase transition were calculated. In general, however, those models fail in their extrapolation 
back to normal matter and lack an adequate description of the groundstate properties of 
nuclear matter, especially they generate too large values for its compressibility. 

On the other hand there has been significant progress in the description of finite nuclei 
using a relativistic meson field (RMF, "Walecka" type) ansatz for modelling the basic degrees 
of freedom of the system |], @]. Those calculations are quite successful in describing nuclear 
properties over a wide range of mass numbers [Q, |]]. 

It seems natural to study model descriptions that combine both of those regimes. There 
has been some effort in the past to follow this path. It turned out that one can get a 
sensible behavior of highly excited matter, a good description of nuclear saturation as well 
as a reasonable description of nuclei and hypernuclei with a single model and a single set of 
parameters. 

In this article we follow and extend this approach. We first determine an improved param- 
eter set || in a fit to spherical nuclei. We then expand our calculations to two-dimensional 
axially symmetric systems. We study several test cases of deformed and superdeformed 
nuclei and perform a general survey of the nuclear quadrupole deformation across the whole 
nuclear chart. We conclude with a critical analysis of necessary improvements to the current 
approach. 
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II. MODEL DESCRIPTION 



In our calculation we use a hadronic model based on a chiral SU(3) ansatz in a non-linear 
realization of chiral symmetry. The basic hadron fields are the SU(3) multiplets for the 
baryons and mesons, respectively, including the baryon octet B that contains nucleons and 
the A, £ and H hyperons, and the corresponding octets and singlets for scalar, pseudoscalar, 
vector and axialvector mesons. In the non-linear realization the pseudoscalar degrees of 
freedom generate chiral rotations of the baryons: 

B = u j ^ L u =u^ R u j . (1) 

with u = exp(iT a \ a /a ). ipi and x/jr are the left- /right handed multiplets in the linear repre- 
sentation of the theory. The 3x3 matrix u contains the pseudoscalar fields 7r a 0] 
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The values <To and Co are the vacuum expectation values of the scalar meson fields as discussed 
below. The renormalization factors containing the ratio w = v^Co/^o are included to obtain 
the canonical form of the kinetic energy terms for the pseudoscalar mesons 0. Note that 
after the transformation (1) the resulting baryons B transform vectorially. 

The baryons couple to the scalar and vector mesons (3x3 matrices W) according to the 
general SU(3) scheme 

C BW = -V2gf (oc w [BOBW] F + (1 - a w )[BOBW] D ) - gf^=Tr{BOB) TrW , (3) 

with [BOBW] F := Ti(BOWB - BOBW) and [BOBW] D := Ti(BOWB + BOBW) - 
^Tt(BOB)TtW and O = {1,7 M } . For the parametrization discussed below we get a ratio 
of the D (anticommutator) and F-type coupling (commutator) of about D/F ~ 0.34. In 
the case of the vector mesons in accordance with vector meson dominance arguments we 
restrict ourselves to pure F coupling. This means that in contrast to the scalar case the 
neutral vector meson made of strange quarks, the (ft ~ (ss) does not couple to nucleons and 
therefore does not enter the equations for nuclei without hyperons. It would be, however, 
phenomenologically consistent to add a small D coupling admixture to the baryon-vector 



meson coupling. The effect of such an admixture should be studied in future investigations. 
The self-interactions of the scalar mesons are taken into account up to 4th order in the fields. 
In addition the interaction terms with the glueball field \ generate a logarithmic potential. 
The field \ mimics the breaking of the QCD scale invariance due to the gluon condensate 
X ~< F^ u F a,flu >. The specific structure of the interaction terms are a direct result of this 
analogy (for a more extended discussion, see [0). 

For the study of nuclear matter and static nuclei in a mean field approximation the 
following fields have to be taken into account: The proton and neutron, the scalar meson 
fields a ~ (uu + dd), ( ~ (ss), the isovector 5 ~ (uu — dd), the scalar glueball field x, and 
the time components of the isoscalar and isovector vector mesons Uq, p[} and the Coulomb 
field A . Restricting ourselves to those degrees of freedom the general SU(3) Lagrangian 
reduces to the following structure. The interaction of the baryons and mesons (and the 
photon) reads 

&nt = ~ Bi fl^oTo + 9i P T3Polo + \e(l + t 3 )v4 7o + m* Bi (4) 

where now B is reduced to the isospinor (£). The various coupling constants of mesons and 
baryons result from the SU(3) structure (|3|). The vector meson mass and self-interaction 
terms are |J 

£ vec = ~\k^ (m> 2 + my) + gi {4 + 6c 2 (p°) 2 + (p°) 4 ) (5) 
and the interaction of the scalar fields follows as 

= -^oxV + C 2 + + h(a" + C 2 + 5 2 ) 2 + h ( £ + £ + 3^ 2 + A 



, 4 1 4, X , 4i 

+hx° C-hx - -.X In— + ex In — (6) 

4 Xo CT oCo 



In addition, the terms 
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£esb = -(JJ [xa + yC] (7) 

break the chiral SU(3) symmetry explicitly with the coupling strengths x = m^f n and 
y = \/2m'^fK — -^ r rnlf 7r , where f w and fx are the pion and kaon decay constant, respectively. 
The baryons attain a dynamically generated mass due to their interaction with the scalar 
fields. This mass generally changes in the nuclear medium. The baryon masses are given by 

m* = g ia o- + g l5 5 + g K ( . (8) 



TABLE I: Set of spherical nuclei used to fit the model parameters of the parameter set Xm- The 
deviations of binding energy, charge radius, and LS splitting values (averaged over proton,neutron) 
are shown. 





6e(%) 


5r ch (%) 




6e(%) 


5r ch { 


48 Ca 


0.64 


-1.38 


124 Sn 


0.24 


-0.23 


56 Ni 


-0.61 


-0.53 


132 Sn 


-0.07 




58 M 


-0.57 


-0.47 


136 Sn 


-0.17 




80 Zr 


0.81 




136 Xe 


0.22 




90 Zr 


0.30 


-0.50 


144 Sm 


0.31 




84 Se 


0.40 




190p b 


-0.06 


-0.96 


88 Sr 


0.29 


-0.64 


202p b 


-0.18 


0.56 


96p d 


0.12 




208p b 


-0.11 


0.46 


ioo Sn 


0.02 




214p b 


-0.17 


0.61 


112 Sn 


0.20 


-0.02 


210p o 


0.02 


1.04 


120 Sn 


0.33 




214 Ra 


0.29 





The asterisk indicates that the masses in the medium shift with the changing isoscalar scalar 
fields a and ( and the isovector scalar field 5. Due to the general SU(3)-symmetric structure 
the S meson automatically enters the field equations for the nuclei including contributions 
to the nonlinear terms of the scalar interaction. Note that in || for simplicity the S meson 
was neglected although this was not quite consistent with the general model approach as 
outlined here. 

III. CALCULATION OF NUCLEI 

We follow the general approach outlined in ||. In contrast to our previous parameter 
search we used a more extended set of nuclei to fit the model parameters as shown in 
Table |. We only used the binding energies of the nuclei for the x 2 fit- The charge radii, 
which generally show reasonable behavior, and the LS splitting (see Table ||) were obtained 
without including them in the fit. Comparing the errors with similar studies in relativistic 
mean field theories || shows that this approach can obtain the same fit quality as well- 
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TABLE II: LS splitting for parameter set Xm- The deviations from experiment are quoted. In 
brackets the corresponding single particle states are stated. 

S(Ls) p (%) 5(LS) n (%) 

16 -13.7 (lp) -9.8 (lp) 

40 Ca -24.1 (Id) -12.4 (Id) 

132 Sn -3.3 (2d) -12.4 (2d) 

208 Pb 0.8 (2d) -29.7 (3p) 



established RMF parameter sets. The fit parameters for the fit Xm is shown in table |T1 



Using the fitted parameters we check that the description of vacuum masses and nuclear 



matter values are still in good agreement with experiment. The result is listed in Table |V. 
Hadron masses changed only slightly. The nuclear compressibility k has a value of 

K = 9p 2 d2 ^( A) =215.33 MeV . (9) 
op 2 

Note that the model produces a satisfactory value for the asymmetry energy defined as 

E asym = \pl S {E,A \,2 = 31-9 MeV , (10) 

1 [d{Pp - Pn)\ 

which is more in accordance with phenomenological values as compared to the usually too 
large values for E asym in RMF calculations ||. 

Here the nonlinear terms in the vector self-interaction (g) and of the scalar S meson are 
crucial for a low value of E asym . In Fig. [I] one sees the influence of the S meson on the 
asymmetry energy. One can observe a softer energy dependence of the system on changing 
the isospin . The total effect is a lowering of the asymmetry energy from ~ 34.5 to 31.9 MeV. 
The new parameter fit and the addition of the 5 field also has some effect on neutron star 
properties due to the general softening of the nuclear equation of state. Compared to an 
earlier calculation [H| the maximum mass of a (non-rotating) neutron star in this model, 
determined by integrating the Tolman-Oppenheimer-Volkov equations as explained in detail 
in flOf , is reduced to M ~ 1.46 M so i ar compared to M ~ 1.64 M so i ar without the inclusion 
of the S meson. 

We check the general applicability of our model by performing a full set of calculations 
of 921 even-even nuclei (first in spherical approximation). The upper panels of Fig. shows 
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FIG. 1: Binding energy of nuclear matter as function of density. Right figure: symmetric nuclear 
matter. Left figure: energy for varying isospin density p p — p n at fixed nuclear matter satura- 
tion density p = .153 GeV/fm 3 . Results with (full line) and without (dotted curve) 5 meson 
contributions are shown. 

the result of a calculation of all even-even nuclei, listed in the Audi-Wapstra-compilation 
jnj (augmented by possible superheavy nuclei up to a charge Z = 136, see Figure below). 
The two-proton and two-neutron gaps are displayed, defined as 

6 2p (N, Z) = E{Z + 2, N) — 2E(Z, N) + E{Z - 2, N) 

6 2n (N, Z) = E(Z, N + 2) - 2E(Z, N) + E(Z, N-2) (11) 

In the plots the average value 52 P (Z) of 52 P (N, Z) for the various isotopes is shown: 

S 2p (Z)^-J25 2p (N,Z) (12) 

71 N 

where n is the number of (even-even) isotopes. The same procedure is performed for the 
neutron gap: 

5 2n (N) = -Y,S2n(N,Z) (13) 

n z 
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FIG. 2: Two-nucleon gaps 52 P (left panels) and #2 n (right panels) in units of MeV for the Xm pa- 
rameter set. The upper panels show the results of a calculation assuming spherical symmetry. The 
lower panels display the corresponding results for a 2-dimensional axially symmetric calculation. 

One can clearly observe the major shell and subshell closures. Looking beyond the tab- 
ulated nuclei in the range of the superheavies results for a spherical calculation is shown in 
Fig. |3|. There is some signal for a shell gap that can be seen in the Id calculation of 52 P for 
Z=120. There is no indication of a gap at the canonical value of Z=114. 

Looking at (upper right panel of Fig. one can observe a weak signal for a shell 
closure at N=172 and N=184, which is in accordance with most relativistic calculations (see 
||12||). Note the strongly suppressed size of the gap compared to normal nuclei as seen in 
Fig. |. 

Let us turn to first 2-dimensional calculations performed in this model. Here the field 
equations are solved on a space-time grid adopting axial symmetry of the system. Part of the 
computer code is an adaptation of the RMF code by Rutz et al.H. As pairing interaction 
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FIG. 3: Two-nucleon gaps 52 P (left panels) and b<m (right panels) in units of MeV for the Xm 
parameter set. Results for hypothetical superheavy nuclei are shown. The upper panels show the 
results of a calculation assuming spherical symmetry. The lower panels display the corresponding 
results for a 2-dimensional axially symmetric calculation. 

we use a zero-range pairing interaction in the same way as it was already applied to RMF 
calculations [IB]. The fitted pairing forces are quoted in Table [TV| The resulting two-nucleon 
gaps in an axially symmetric 2-dimensional calculation are shown in the lower panels of 
Fig. |2|. The general behavior is similar to the spherical results. The proton subshells are 
suppressed, the major shell closures are clearly visible. The situation for neutrons is similar, 
there is still a signal of a subshell at iV = 40. If we look at the superheavy elements in 
the 2-dimensional calculation (lower panels in Fig. |3| the signal is completely washed out, 
and at least the gap energy cannot serve as a guideline for identifying shell closures. This 
originates from the strong deformation of the nuclei involved in calculating equations (11). 
In the equation nuclei with different deformations enter making it impossible to directly infer 
conclusions on the shell structure of the nuclei. This is similar to a possible shell quenching 
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effect at Z=82 that is largely artificially generated by comparing deformed nuclei above and 
below the shell closure as discussed in |12| . 



Comparing the 2d calculations with other theoretical calculations we study the chain of 
magnesium isotopes. Here we perform a calculation with a constraint on the quadrupole de- 
formation of the nucleus (leaving all other multipoles unconstrained). We vary the constraint 
and calculate the energy of the isotopes as function of their deformation. The resulting fig- 
ure is shown in ||. These results that show the occurrence of oblate and prolate minima for 
various isotopes is in agreement with other relativistic and non-relativistic mean-field calcu- 
lations [14], |15, [16], [17| . There is no strongly deformed groundstate for 32 Mg as experimental 



measurements of B(E2) values suggest |[18|| . One can, however, observe a shoulder in the 
energy at deformations /?2 ~ 0.1 to 0.2, which is due to an intruding f7/2 neutron state. 

The overall quality of the model for describing finite nuclei can be estimated from the 
error in the binding energy averaged over all nuclei. The comparison has been done using the 
Audi-Wapstra experimental data table. As a total average for nuclei with charge Z= 8 and 
higher we obtain a deviation e# ~ 0.35% from the the experimental binding energy values. 
Note, however, that most deviations arise from light nuclei A <= 40 where more subtle 
effects affect the binding energy that cannot accurately be reproduced in a mean field-type 
description as applied here. Considering only heavier nuclei we get the following numbers: 

e B m (A >= 50) = 0.21%, e B m {A >= 100) = 0.14% (14) 

We can compare this result with an analogous calculation in the RMF model (using the NL3 
parameter set with zero-range pairing force [|13|j ) 

e N B L \A >= 50) = 0.25%, e^ L3 (A >= 100) = 0.16% (15) 

We conclude that the quality of the nuclear fit within our model is slightly better quality 
compared to standard relativistic nuclear models. Comparing the results with macroscopic- 
microscopic calculations of Moller and Nix, using their data compilation [ft"9|l , we read off 
values of e MN (A >= 50) = 0.05% and e MN (A >= 100) = 0.03%. This shows that in self- 
consistent models one is still off by a factor of about 4 to 5 from those numbers (Note that in 
contrast to [T^] the numbers from \ m and NL3 arise from purely axially and time-symmetric 
calculations). 

We calculated the nuclear quadrupole deformations over the whole range of known nuclei 
assuming axial symmetry. The calculation was done in the following way. We initialized each 
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FIG. 4: Energy (relative to the spherical calculation) of Mg isotopes as function of quadrupole 
deformation 

nucleus starting with 3 configurations - prolate, oblate, and spherically symmetric and let 
the calculation converge. We then took the deformation of the energetically lowest solution 
as groundstate deformation of the nucleus. To a large extent this procedure reduces, but 
it does not completely eliminate, the chance of missing out on the true minimum because 
of shape isomers . There is still the possibility that a energetically lower state has been 
missed. One can improve on this procedure by doing a constrained calculation for every 
nucleus, varying the quadrupole deformation and picking out the state with lowest energy. 
We will look at some cases below. Obviously, this method needs much more computer time. 
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FIG. 5: Groundstate quadrupole deformation |^| hi our model using the Xm parameters in a 
2-dimensional calculation assuming axially symmetry. White squares denote spherical shapes, the 
darkest shade of gray represents deformations with \(3%\ > 0.4. 

Here, the total time used for generating Figure |5] was about 2500 CPU hours on a 500MHz 
Pentium machine. 

There are a number of areas where shape isomers exist. Especially in the Z ~ 70 range 
and large neutron number the calculated energies of the oblate and prolate minima are often 
within a few hundred keV, and, given the approximations involved in this calculation, the 
true groundstate deformation cannot really be determined (in a more realistic calculation us- 
ing configuration mixing both minima will mix). It is possible that both states are connected 
via 7 deformation. In future 3-dimensional calculations this point should be studied more 
closely. As an example the Fig. shows the resulting deformation dependence of the energy 
of very neutron-rich Dysprosium and Erbium isotopes. You see quite deformed prolate and 
oblate minima with f3 2 ~ 0.45 and —0.28, respectively, that are essentially degenerate. In 
principle, there are two values for the deformation to be displayed in Fig. || 

In a comparison with recent experimental results of the deformation of various sulfur and 
argon isotopes |20| we calculated the groundstate deformation of 38_46 S and 42,44 Ar. The 
results can be seen in Figure [7]. The experiments determine the absolute value for 02, the 
calculations show prolate groundstates for sulfur and oblate ones for argon. Our numbers 
agree quite well with experiment, some predictions for the sofar undetermined deformation 
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FIG. 6: Binding energies of Dy and Er isotopes, relative to the spherical state, as function of defor- 
mation. Two nearly degenerate minima can be observed. The horizontal lines show experimental 



(or extrapolated) values for the binding energy [11| 



of very neutron rich sulfur isotopes are shown, too. Going to proton-rich side we compare 
our results of a calculation of axially deformed N=Z nucleus 68 Se with relatively recent 
experimental data that show a strongly oblate groundstate (fa ~ —0.3) together with an 
excited prolate state [2~I|. The result of a quadrupole-constrained calculation is shown in 
Fig. [| In accordance with the experiments distinct strongly deformed oblate and prolate 
minima are found, where the oblate state is more strongly bound by about 2 MeV. As a 
further test of calculations of deformed nuclei we look at Nobelium isotopes. The calculation 
of the deformation of 252 No, 254 No, and 256 No are shown in Figure ^. The behavior of the 
different isotopes is very similar, showing a prolate minimum at around ~ 0.29. This 
compares very well with measurements |22| of the groundstate deformation of 252 No and 
254 No with resulting numbers of (3 2 = 0.31 ± 0.02 for A=252 and fa = 0.32 ± 0.02 for A 
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FIG. 7: Quadrupole deformation |/?2| for sulfur (circles) and argon (squares) isotopes with mass 
number A. The full symbols denote the theoretical results, the open symbols are experimental data 



TABLE III: Parameter set Xm, ah numbers are quoted for the appropriate powers in MeV. G p and 
G n are the strengths of the zero-range pairing force for protons and neutrons, respectively. 



9Na 


-10.569 


9N( 


0.467 


9Nw 


13.3265 


9Np 


5.48851 


k 


2.47584 


fci 


1.35436 


k 2 


-4.93719 




-2.77257 




-0.233064 


94 


79.9083 


e 


0.02134 


Xo 


409.76956 


G p 


-4.47- 10~ 5 


G n 


-4.34 • 10~ 5 


fK 


122 







= 254, respectively. Finally, we investigated various super deformed heavy nuclei. In all 
cases where a super deformed state is known we could either see a distinct superdeformed 
minimum or a shoulder indicating the precursor of a superdeformed state for higher spins. 
As an example Figure [T^ shows the result for the mercury isotopes 192 Hg and 194 Hg, which 



has originally been predicted to have superdeformed bands in p3| (a study analogous to our 
discussion using in the RMF approach has be performed in [24]). In addition to the oblate 
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FIG. 8: Binding energy of the nucleus 68 Se (relative to the spherical solution) as function of the 
quadrupole constraint /?2- 

TABLE IV: Hadron masses (in MeV) for parameter set Xm- In brackets the experimental isospin- 
averaged values are shown p5[ . Note that in the particle data book a, £, <5 correspond to /o(400 — 
1200), / (980), and o (980), respectively. 

m N 939.2 (938.9) m A 1115.7 (1115.6) m 2 1196.0 (1193.1) m~ 1331.3 (1318.1) 
m n 138.6 (138.0) m K 497.3 (495.6) m v 574.0 (547.5) 897.3 (957.8) 

m a 466.5 (400 - 1200) m c 1024.5 (980) m 5 973.26 (985) 1020.0 (1019.4) 

m u 780.6 (782.0) m p 761.1 (768.1) m x 960.0 

groundstate we see a rather complex structure of the energy surface and clear evidence 
for super deformation in both cases. An experimental measurement and extrapolation to 
the groundstate of 194 Hg gave an estimate of an excitation energy E ~ 6.02 MeV |26] of 
the superdeformed state with respect to the oblate minimum. In our calculation we get 
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FIG. 9: Binding energies of different Nobelium isotopes as function of the quadrupole constraint 
/?2- The horizontal lines indicate the measured groundstate energy 

E = 4.3 MeV, which is a little bit low but in reasonable agreement with experiment. In the 
case of 192 Hg we obtain E — 3.1 MeV. 

IV. CONCLUSIONS 

We investigated a hadronic model based on chirally symmetric interactions within an 
flavor-SU(3) approach and calculated nuclear properties within this approach. First we 
introduced a new improved set of model parameters generated by a fit of a set of spherical 
nuclei. For the first time also the isospin-triplet 5 meson was taken into account, which 
naturally occurs within the SU(3) scheme in its coupling to the baryons as well as in the 
nonlinear meson-meson interactions. The quality of the fit is as good as the results of other 
relativistic mean field models of nuclear structure. We investigated deformed nuclei and 

16 



-1510 



*£ -1516 

m 



-1520 



192 




■ ■■ 1 1 1 1 ■■■■■■■■■ 1 1 1 1 1 ■■■■■■■■■ 1 1 1 1 ■■■■■■■■■ i 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 




1.0 1.2 



FIG. 10: Binding energies of 192 Hg and 194 Hg as function of fa- Both isotopes exhibit an oblate 
groundstate and super deformed structures for ~ 0.7. Experimental binding energies are indi- 
cated by horizontal lines. 



found reasonable results for nuclear quadrupole deformations. A number of improvements 
should be addressed in future calculations. A projection onto good particle numbers as well 
as, in case of the deformed calculations, a projection onto good angular momentum should 
be performed. In different regimes of the nuclear landscape shape isomers are found with 
minima that are very close in their energies. Especially for those cases three-dimensional 
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studies and a calculation including configuration mixing of the various states have to be 
done. All these points are currently under investigation [27| . 
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